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Abstract. A pair of linearly independent asymptotic solutions are con- 
structed for the second-order linear difference equation 

Pn+l{x) - {AnX + B„)Pn{x) + Pn-l{x) = 0, 

where A^ and i3„ have asymptotic expansions of the form 



s=0 s=0 



with 7^ and ao 7^ being real numbers, and /3o = ±2. Our result 
hold uniformly for the scaled variable t in an infinite interval containing 
the transition point t\ = 0, where t = (n + to)~^x and tq is a small shift. 
In particular, it is shown how the Bessel functions Ji, and Yy get involved 
in the uniform asymptotic expansions of the solutions to the above three- 
term recurrence relation. As an illustration of the main result, we derive 
a uniform asymptotic expansion for the orthogonal polynomials associated 
with the Laguerre-type weight a;" exp(— 9^2;™), a; > 0, where m is a positive 
integer, a > — 1 and (7,„ > 0. 
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1. Introduction 

Orthogonal polynomials play an important role in many branches of math- 
ematical physics, for instance, quantum mechanics, scattering theory and sta- 
tistical mechanics. A major topic in orthogonal polynomials is the study of 
their asymptotic behavior as the degree grows to infinity. Since the classical 
orthogonal polynomials (Hermite, Laguarre and Jacobi) all satisfy a second- 
order linear differential equation, their asymptotic behavior can be obtained 
from the WKB approximation or the turning point theory [15]. For discrete 
orthogonal polynomials (e.g., Charlier, Meixner and Krawtchouk), one can use 
their generating function to obtain a Cauchy integral representation and then 
apply the steepest descent method or its extensions [151 [23]. However, there are 
orthogonal polynomials that neither satisfy any differential equation nor have 
integral representations. A powerful method, known as the nonlinear steepest 
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descent method for Riemann-Hilbert problems, has recently been developed 
that can be applied to such polynomials. Papers that deserve special mention 
include Deift et al. [9], Bleher and Its [6], Kriecherbauer and McLaughlin [H], 
Baik et al. [2], Ou and Wong [17] and Zhou et al. [27]. Despite the great 
success achieved by this powerful method, the Riemann-Hilbert analysis de- 
pends heavily on the analyticity of the weight functions, and the argument and 
results obtained by the method often appear in a very complicated manner. 

In our view, a more natural approach to derive asymptotic expansions for 
orthogonal polynomials is to develop an asymptotic theory for linear second- 
order difference equations, just in the same way as Langer, Cherry, Olver had 
done for linear second-order differential equations; see the definitive book by 
Olver [15]. Our view is based on the fact that any sequence of orthogonal 
polynomials satisfies a three-term recurrence relation of the form 

(1.1) Pn+l{x) = {anX + hn)Pn{x) - CnPn-l{x), 72 = 1,2,..., 

where a^, 6„ and constants; see [TBI P-43]. If ^ IS cL fixed number, 

then this recurrence relation is equivalent to a second-order linear difference 
equation of the form 

(1.2) y{n + 2) + n^'a(ra)y(n + 1) + n%{n)y{n) = 0, 

where p and q are integers. We assume that the coefficient functions a{n) and 
b{n) have asymptotic expansions 

oo oo , 

(1.3) «(")~E^. 

s=0 s=0 

where ao 7^ and 60 7^ 0. When p = q = 0, asymptotic solutions to this 
equation are classified by the roots of the characteristic equation 

(1.4) + aop + bo = 0. 

If pi 7^ p2, i-G-, 7^ 4605 then Birkhoff [3] showed that (11. 2p has two linearly 
independent solutions, both of the form 

00 

(1.5) ?/(n) -p'^n"^^, n^oo, 
where 

aip + bi aip + bi 

(1-oj a = = — , 

2p^ + aop 2bo + aop 

and the coefficients Cg can be determined recursively. This construction fails 
when and only when pi = p2, i.e., when Og = 46o. Motivated by the termi- 
nologies in differential equations [151 P-230], we shall call series of the form 
(11. 5p normal series. If pi = p2, but 2bi 7^ Ooai, then Adams [1] constructed 
two linearly independent power series solutions, which we will call subnormal 
series. 



LINEAR DIFFERENCE EQUATIONS WITH A TRANSITION POINT AT THE ORIGIN 3 



Adams also studied the case when 2bi = aoOi, but his analysis was incom- 
plete, as pointed out by Birkhoff A powerful asymptotic theory for dif- 
ference equations was later given by Birkhoff Birkhoff and Trjitzinsky [5]. 
However, their analysis has been considered too complicated and even impen- 
etrable. A more accessible approach to the results mentioned above has been 
given in more recent years by Wong and Li [211 [25] . In their first paper [2l] , 
they presented recursive formulas for the coefficients in the power series so- 
lutions when p = q = in fll.2p . obtained both norm and subnormal series 
solutions. Their second paper [25j dealt with the more general case, namely, 
equation f ll.2p with the exponents p and q being not both zero. Numerically 
computable error bounds for these asymptotic series can be found in our recent 
paper [7]. 

Returning to the three-term recurrence relation (11.11) . we note that the re- 
sults mentioned thus for apply only when a; is a fixed number. When a; is a 
variable and allowed to vary, the roots of the characteristic equation associated 
with (11.11) may coalesce, and normal series solutions may collapse to become 
subnormal series. Not much work has been done in this area until just recently. 
In a series of papers [2niEIll22j, Wang and Wong have derived asymptotic ex- 
pansions for the solutions to (11.11) . which hold "uniformly" for x in infinite 
intervals. They first define a sequence {K^} recursively by K^+i/Kn = Cn, 
with Kq and Ki depending on the particular sequence of polynomials, where 
Cn is one of the coefficients in equation (II. ip . Then they put An = anKn/ Kn+i, 
Bn = bnKn/Kn+1 and Pn{x) = pn{x)/Kn-, SO that (II. ip becomes 

(1.7) Pn+l{x) - {AnX + Bn)Pn{x) + P„_i(x) = 0. 

The coefficients An and Bn are assumed to have asymptotic expansions of the 
form 

oo oo „ 

(1-8) -^"-""E^. 

s=0 s=0 

where ^ is a real number and ao 7^ 0. If Tq is a constant and N := n + tq, then 
the expansions in (II. 8p can of course be recasted in the form 

00 , 00 „/ 

(1-9) A.~iv-"E^. «»~E^- 

s=0 s=0 

In ([LID, we now set x := NH and P„ = p". Substituting ([L9]) into ([L7|) and 
letting n — > 00 (and hence — > 00), we are led to the characteristic equation 

(1.10) p"" -{a'ot + /3'o)p + 1 = 0. 

(Note that «[, = ao and /3q = (3o). The roots of this equation coincide when 
t = ti, i = 1, 2, where 

(1.11) a'oti + (3'o = +2, a'ot2 + /3'o = -2. 
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The values of ti and t2 play an important role in the asymptotic theory of 
the three-term recurrence relation fll.?!) . and they correspond to the transition 
points {i.e., turning points or singularities) occurring in differential equations 
[T5] . For this reason, Wang and Wong [21] also called them transition points 
of difference equations. Furthermore, they indicated that in terms of the ex- 
ponent 6 in (11. 8p and the transition point ti, there are three cases to be con- 
sidered; namely, (i) 6' 7^ and ti 7^ 0, (ii) 6^0 and ti = 0, and (iii) ^ = 0. In 
[21], they considered case (i), which turns out to correspond to the turning- 
point problem for second-order linear differential equations [T^ p. 392]. They 
also showed that the asymptotic expansions of the solutions involve the Airy 
functions Ai(-), Bi(-) and their derivatives, and furthermore that the expan- 
sions hold uniformly for t in the infinite interval [0, 00). In [22], they studied 
case (iii), and found that the approximants of the asymptotic solutions are 
Bessel functions or modified Bessel functions. We should mention that, the 
WKB approximations for difference equations have also been studied in recent 
years by Costin and Costin [8], Geronimo, Bruno and Van Assche [13], Geron- 
imo [12], etc. All of these authors dealt with only case (i), i.e., the Airy- type 
expansions. 

The purpose of this paper is to provide a solution to the problem in case (ii), 
which is the only case that have been left out in Wang and Wong's investigation 
of the asymptotic solutions to the three-term recurrence relation (11. 7p . In this 
case, we shall show that the approximants of the asymptotic solutions are also 
Bessel functions or modified Bessel functions, but the order of these Bessel 
functions and the transformation used {i.e., ({t) in (14.91) below) differ from 
the one used in case (iii). As an illustration of the main result, we derive a 
uniform asymptotic expansion for the orthogonal polynomials associated with 
the Laguerre-type weight exp(— g^^^™"), a; > 0, where m is a positive integer, 
a > —1 and > 0. 



2. Motivation leading to the expansion 

Turing point (or transition point) theory for second-order linear difference 
equations was well estabhshed in [2T1 [22] • We simply state the procedure 
below. 

Let To := —ai/{aQ9) and N := n + Tq. Then, 

(2.1) ^„. + B„ = „-.^|. + j;^:=^f!£±A. 

s=0 s=0 s=0 

A slight computation yields 

(2.2) a'o = ao, a[ = 0, (3'^ = (3o, (3[ = (3,, = + /^iTq. 

If Pq = 2 (or —2), it follows from (II. lip that one of the transition points is 
zero. Without loss of generality, we may assume that ti = < ^2 and /3o = 2, 
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cto < 0. (For other cases, see Remark [T] below.) Through out this paper, we 
assume that /3i = 0, so that /3[ = and 



(2.3) 



a'lh + /3( = /3i = 0. 



This assumption was used in the previous papers on transition point theory; 
see [21], (2.7)] and |22l (2.5)]. (In most of the classical cases, /3i = and hence 
f3[ = 0. Also, we note that in [10] Dingle and Morgan have assumed a more 
strict condition that a2s+i = (32s+i = for s = 0, 1, 2 . . . .) 

In what follows we will investigate the uniform asymptotic expansions of 
solutions to (11. 7p near the transition point ti = 0. To this end, we seek a 
formal solution to (11. 7p of the form 



(2.4) 



s=0 



for t in a neighbourhood of ti = 0, where r]{t) is an increasing function with 
?7(0) =0. Since x := N^t is fixed, when we change n to n ± 1, we must in the 
same time change t to t±, where 



1 ± — 

N 



t. 



(2.6) Pn±iix) = J2x. 

s=0 



(2,5) (± 
As a consequence, it follows from (12. 4p that 

{N±iy7] 

For convenience, we also introduce the notations 
(2.7) Q±{0:={N±l)'v 
and 
(2.8) 



1 





-e ' 


[04) 


t 







s=0 



Substituting (E^]) and ([M]) into ([LTD, we get 
(2.9) 



f 1— 1— 'ii 



Noting that t = rj ^[N ^^), Taylor expansion gives 
^ r/"^(^)(0) 1 ^ 



(2.10) 



j=0 



m J r/'(0) 
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Substituting the last equation into (12. 8p yields, 

(2.11) To(0 = = /3o, Ti(O = /3;=0, r2(0 = ^ + /3^. 

By expanding ri(t±) at t, we obtain from (12. 5p and (12. 7p 



In equation (12. 9p . the coefficient of A^~* vanishes for each s > 0; i.e., 

(2.13) Xs[Q+m (i + ^)"' + Xs[Q-m (i - ^)"' - = 0. 

In particular, for s = 0, expanding the ffist two terms at ^ and using (12. 8 p and 
(I2.12P show that Xo(0 satisfies Bessel's equation 

where 6 2. For simplicity, we shall assume < < 2. The analysis for the 
case ^ < or > 2 is very similar. Thus, Xo(0 can be expressed in terms of 
Bessel functions: 



where Ci, C2 are two constants, b = \/ (g_2)^"^°(o) and 



(2.15) 



1 + 4/3^ 



Moreover, each of the subsequent coefficient functions XsiO^ s = 1,2, . . . , in 
(12. 4p satisfies an inhomogeneous Bessel equation. This suggests that instead 
of (12. 4p . we might try the formal series solution 

(2.16) p„(x) = z.(ivc) f: ^ + ^.+i(ivc) f: ^ 

s=0 s=0 

motivated from the differential equation theory, where we have set ({t) = 
[r){t)]2 to simplify the notations. Here (^{t) > for t > and iC(t) < for 
t < 0. In (I2.16p . Zu{-) could be any solution of Bessel's equation 

(2.17) y" + -y'+ (l-'^]y = 0. 



We state the main result of this paper in the following theorem. 
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Theorem 1. Assume that the coefficients An and Bn in the recurrence rela- 
tion ( [j. ?] ) are real, and have asymptotic expansions given in U.8\) with 9 0,2 
and = 2. Let ti = be a transition point defined in U.ll\) . and the function 
({t) be given as in Ili4.9\ ) and Then, for each nonnegative integer p, 

( [i. ?| j has a pair of linearly independent solutions. 



(2.18) 



and 



(2.19) 



4e 



4 - ia',t + 



UNO 



MO 



Qn{Nh) =m 



4C 



4 - {a'^t + 



'\2 



s=0 



p 



+ W, 



(NO J2^ + 6^n 



s=0 



where W^{x) := Y^^x) —iJ^{x), N = u + Tq, Tq 
\2.15\) . The error terms satisfy 



-ai/(ao^) (^''^d V is given in 



(2.20) 

and 
(2.21) 



M 



'-A\MNO\ + \Ui{NO\] 



for —oo<t<t2 — (y, Mp being a positive constant and a being an arbitrary 
positive constant. The coefficients v4o(C) = 1 o^nd Bq{() = and other As{Q 
and -Bs(C) can be determined successively from their predecessors; see ( [^^ . 26] ) 
and [jlTD . 

Remark 1. For the case > and Po = —2, the two transition points 
ti and t2 again satisfy ti = < ^2- Put Vn{x) := {—l)^Pn{x) . Theorem\J\ 
then applies to Vn{x). If the two transition points satisfy t2 < = ti, we set 
X := —N^t, instead of x := N^t. Theorem\l\ again holds with Pn{N^t) replaced 
by Pn{-NH) (m the case ao > 0, f3o = 2) or (-l)"P„(-A^^t) (in the case 
«o < 0, /3o = -2). 



3. A PRELIMINARY LEMMA 



As we have noted in above sections, x := N t is fixed in equation (11.71) . 
Hence, as N being replaced by A^±l, the two functions Z^{NC.) and Zi,^i{N() 
in l^lEf automaticaUy change to [{N ± l)C{t±)] and Z^+i [{N ± l)C{t±)], 
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where t± are given in f l2.5p . The following lemma plays a crucial role in the 
derivation of the formal series solution in fl2.16p . 

Lemma 3.1. Let Z^{x) be any solution of Bessel's equation 

(3.1) y" + -y'+(i-'^)=o. 



Then, we have 

(3.2) {{N ± l)C(t±)} = Z,(iVC)G±(c, ^) + Z,+i(iVC)if±(c, ^) 



and 



(3.3) Z,^, {{N ± l)C(t±)} = Z,{NOL± (C, ^) + ^.+i(iVC)K± (C, ^) : 



where 



(3,4) G,(c.l)^f;(±ir^, ,f,(c,l)^f;(±i)-^ 

s=0 s=0 



and 



.i<.(Q 



^■5) i.(4)~E(±ir'^. '^■44)~E(±i) . 

s=0 s=0 

t/ie expansions being uniformly valid with respect to bounded t. 
Proof. Let 

(3.6) w{N, C, u)=(l + ^y Z, [{N + u)C] . 

A straightforward calculation gives 

(") l:-i^(^-^)"'-'-(-^)'c-- 

and 

where we have made use of the difference-differential relations of Bessel func- 
tions 

Ki^) = -Zu{x) - Z^^i{x) 

X 

and ^ 

Z'y+ii.^) = Zv{x) Z^+i{x)] 

X 

see [m §10.6]. The Taylor expansion of w{N,(,u) at m = is 

(3.9) wiN, C, u) = wiN, C, 0) + n— (AT, (, 0) + -^(AT, (, 0) + ■ ■ • . 
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In view of (13. 6p - ( 13781) . the above series can be rearranged as 



U\2 



(3.10) + Z4{N + u)C] = Z,{N0G{N;C,u) + Z,+^{N0H{N;Cu) 

where 
(3.11) 



G(iV;C,0) = l, H{N-C,0) = 0, 



(3.12) 



dH 



i and :^(iV;C,0) = -C. 



Differentiating (I3.10p with respect to u twice yields 



(3.13) 

and 
(3.14) 



dv? 



2 1 



(A^ + n)2 



dG 
du 



u=0 



2 1 



{N + u) 



1 H, H\u=o = 0, — — 



dH 

du 



+ V 



<■ 



M=0 



We try formal solutions to (13.131) and (I3.14p of the following forms 

oo 

(3.15) G{N-CM-Y. 



and H{N;C,u) ^ 



s=0 s=0 

where the coefficients Gs can be determined recursively by the equations 



(3.16) 
(3.17) 

(3.18) 

and 
(3.19) 



-C^Go, Go(0,C) = l, ^(0,0 = 0, 



du 



||i^-C^g. 5.(0.0^0. ^(0.0 4 + ^ 



3'G, 

du^ 



-eG,+ U-^^ j^{3~i){-ur'G, 



~ dG 

a(0,C) = 0, —^(0,0 = 0, s>2. 



-3-> 



s > 2, 



du 



The solutions of (I3.16p - (l3.19p are 



(3.20) 

and 

(3.21) Gs 



Gq = cos(Cm), Gi 



+ V 



sin(CM) 



^2 _ 1 ru , s 



Y.{-\)K2 - 1)0^"'^-,) sin ((« - 0)C)d(/. 

i=2 
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for s > 2. Similarly, we also have 

(3.22) ^o = -sin(CM), = 
and 

(3.23) = ""—^ \ (5^(-l)^(j - 1)0^-'^.-,) sin {{u - 0)C)ci(/. 

for s > 2. By induction, we can prove that when i^ < 0, 

|a(n,C)| < + l)>rGo(w,C), \tls{uX)\ < iW\ + l)>|1^o(w,C)|, 
and that when C > 0, 

\Gsiux)\ < iw\ + i)>r, \Hsiu,c)\ < iw\ + i)>r- 

Thus, the formal series in fl3.15p are uniformly convergent for any bounded u 
and sufficiently large A^. 

To show that (13.21) follows from (I3.10p . we simply set 

(iV + n)C = (iV±l)C(t±) 

i.e., we choose 

(3.24) u = u±:={N±l)^-N. 
Furthermore, we let 

(3.25) G± (C, ^) = (l + ^) ' GiN; C, u^) 
and 

(3.26) = (i + ^)"'^(iV;C,«±). 

Recall that t± is given in (12. 5p . Expand t± into power series in j^, and then 
expand C(''^±) i^ito Taylor series at t and rearrange them also into power series 
in From (Km . it follows that 

oo 

(3-2T) = E(±l)^"'^' 

s=0 

where 

(3.28) .o^i-''*^, 

^ ' " C(i) 2 Cit) 2 C(^) 

Substituting (13.270 into (I3.15p . we obtain 

(3.29) «,(c,i)^(i,^)-f:to.f:^|^, 

s=0 s=0 
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(3.30) h4c.^)^(i.^)-'i:%^-i:^'''' 

s=0 s=0 



Simple calculation shows that the first two terms in the expansions of fl3.29p 
and fl3.30p are, respectively, given by 

(3.31) Go^(C) = cos(Cwo), G^iO = ±- sin(Cwo) T cos(Cwo) 



and 



(3.32) Hi{C) = Tsin(CMo), HfiC) = -miCcos(Cmo) + y sin(CMo). 

Equation (13. 3p and the expansions in (13. 5p can be established in similar 
manners. Moveover, we have 

(3.33) K^{C) = cos(Cwo), KtiO = ^1±1±^^ sin(C«o) T y cos(Cwo) 
and 

(3.34) L^iC) = ±sin(CMo), L±(C) = miCcos(Cmo) - ^sin(CMo). 
By mathematical induction, one can show that 

(3.35) G;(C) = i-iYGtiO, h;{0 = {-ly-'HtiO 

and 

(3.36) k;{0 = i-iYKnO, l;{0 = {-ly-'LtiC). 

When we consider the case '+', for simplicity we shall drop the subscripts and 
superscripts in all above equations. Thus, (I3.35P and (I3.36P give 

(3.37) G,(ci)=G(c.4), ff,(cl)=±i/(c,4 

and 

(3.38) A-,(ci)=/<C.±^). L.(c.^)=±i(c.4; 

This completes the proof of Lemma 13.11 □ 

4. Formal asymptotic solutions 

Let C^(^) = vit) be an increasing function with C(0) = 0; see the sentence 
following (I2.16p . We try a formal series solution to (II. 7p of the form 

(4.1) P^NH) = Z^iNO £ ^ + ^.+i(iVC) £ 

s=0 s=0 
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where is a solution of Bessel's equation. For convenience, we put 



s=0 

and 



Ml 



s=0 



Ml 



(4.3) 



^ oo 

^'(t,-) := A„x + B„ = ^ 

s=0 



«;iM. 



cf. (12.81) . Using (12.161) and Lemma [3. ![ we have 

P„±i(x) =Z„(iVC) {G(C,±^)'l(C(«±),;jy^ 



(4.4) 



±L(c.±i)fl(cfe),^)} 

+ z..i(ivC){A-(c.4)^(c(*±)'ivTT 
±i/(c.4)A(cfe).-^)}; 



see also (13:371) and (13:38|) . Substituting (il])-(l43D into (JLT]) and matching the 
coefficients of Zy and Z^^+i, we get 



(4.5) 



and 



(4.6) 



1 

iV 

-1 



A 



at-) 



1 



iV+ 1 

1 

' A^- IJ 



-4c 



5 



B 



C(t-), ^ 



A^- IJ 



■*(c-^)-^ 



1 



C(t- 



A^+ 1 
1 

' A^- 1 



5 



C(t+), 



at-) 



N+li 
1 



A^- 1 



('■^)Kc.^)^o. 



Letting A^ — )■ oo in (14.51) and (14. 6p . we obtain 



(4.7) 



Go(C) = i^o(C) 



which combined with (I3.28P and (I3.3ip yields 



(4. 



C(t) - OtC{t) = iarccos 
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where the arccos function is analytically continued to C \ {(0, — ioo) U (^2, oo)} 
so that for t < 0, 



arccos 



-log 



+ iWl 



where t2 is given in fll.lip . Recall that < and /3q = 2; see the lines between 
(12.21) and fl2.3p . It is straightforward to verify by direct differentiation that the 
function 



(4.9) 



±C(t) = arccos 



a 74-(a^0 + /3^)2 



is a solution to (14. 8 p for t < ^2 — cr, where '+' sign is taken for 6 <2 and '— ' 
sign is taken for 6 > 2. The choice of the lower limit of integration is just for 
the purpose of convergence of the integral. For instance, we can chose a = if 
< or > 2 and 

r - oo, t < 

(^■i"' "=|*.. a<t<u, ifo<«<2. 

With this choice, it can also be shown that ^(0) = 0. 

Equating coefficients of like powers of in (14. 5 p and (14. 6 p to zero, we obtain 



(4.11) 



s<p, seven j+m<s ^ 

o <'" frl o i-l T _L rvi <^ c ^ 



p + s 
s — i — m 



-p + s 
s — i — m 



s<p,s odd i+m<s 

iEK* + ^^) '»,-,(C) = o 



E 

1=0 

E"' 

1=0 



1 ll,m 



i ll,m 



s<p 



and 



(4.12) 



E E '<'io( 

\p, s even i+m<s ^ 

E E^^<(c)( 

'ti c r\rlrl T -I- rr) <' c ^ 



p + s 
s — i — m 



-p + s 
s — i — m 



s<p,s odd i+m<s 

iEK' + ^^)-8p-(C) = 0, 



E 

1=0 

E"' 

1=0 



D'Ap-s{C) ,i 
1\ ^^'^ 



s<p 



where denotes the j-th derivative with respect to t, i.e., 



(4.13) 



D^MO = —Mat)). ^ = 0,1,2, 
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and 7/,m is defined by 



(4.14) Yl 
For convenience, we define 



1 + 



(4.15) 



2<s<p 

-E 



2<s<p 
s even 



2<s<p 
s odd 



E 

i+m<s 

E 

i+m<s 



-p + S 
s — i — m 

-p + s 
s — i — m 



1=0 



1=0 



and 



2<s<p 



(4.16) 



-E 

2<s<p 
s even 

-E 



EC:r_:)-.(c)E 

i+m<s ^ ' 1=0 



T, lUmt 



E 



2<s<p \_i-\-m<s 

s odd 



-p + s 
s — i — m 



1=0 



n 



for p > I. Clearly, /o(t) = g^it) = 0. N o te tha t 70,0 = 1, 7o,m = 0, m = 
1, 2, ... , and 71,0 = 0, 71,1 = -6. By (I33B - (I33D and (031), equations fliTT]) 
and (I4.12P can be rewritten as 

(4.17) [(1 -p)Lo(C) + Li(C)]i?p-i(C) - etL,{C)j^B,_,{C) = /p_i(t) 
and 

(4.18) [(1 - p)Ho{0 + i/i(C)]A,_i(C) - etHoiOj^A^^iiO = g,^,{t), 

where we have made use of the assumption that ai = /3[ = 0; cf. (12. 2p and 
fra . It follows from fl33TD . fl3:32|) and (glD that 



(4.19) 



HoiO = - 



4-(a[,t + /3^)^ dHoiO _ a'o GpiQ 
A ' dt 2 HoiC)' 



Differentiating (14. 8p . we get 



(4.20) 



:i - 9)Cit) - etC{t) 
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Combining (Km . ( ]332|) . (Km . ( ICTll and (Km . we obtain 



(4.21) 

and 
(4.22) 



l-dt- 



c 



c 



^o(C) 



^o(C). 



Substituting flOT]) and (KT3f into fHTTD and fITO]) and noting that Lo(C) 
—Ho{(), we have 



(4.23) 

where 
(4.24) 



|[tVA(t)i?,_i] 



A(t) :=t5 



t^-^A(t) 
^^o(C) 
t"i'-^A(t) 



c 



For p = I, since fo = go = 0, AAo(C) and ASo(C) si's constants, we set 
(4.25) AAo(C) = 1, ABo = 0. 

For p > 1 and < < 2 



(4.26) 
and 
(4.27) 



t"^A(t)i?p_i(C) 



t^A(t)A,_i(C) 



p-i 



^-iA(s) 

eHo(cRy 



fp-i{s)ds, t<t2-o 



' sV-iA(s) 



gp-i{s)ds, t<t2-a. 



^^o(C(^)) 

Therefore, for each p > 1, Ap(t) and -Bp(t) can be determined successively from 
their predecessors ^0(^)5 -Bo(t), Ap_i(t) and Bp_i(t). If, in addition, we 
have the assumption a2s+i = /^L+i = for s = 0, 1, . . . , we can choose the 
solutions such that A2S+1 = = 0. 

5. Bounds for coefficients 

Lemma 5.1. Let ({t) be given as in ^.9^ and ^.KJj ), and let G{(,j^), 
H{(, j^), L{(, j^) and K(^(, j^) be given as in Lemma l3J[ Then there exists a 
constant Cg independent of t such that 

(5.1) 1^(01 < C,i\t\ + 1), \HsiO\ < Csi\t\ + 1), 

(5.2) \KsiC)\ < Csi\t\ + 1), |L,(C)| < Csi\t\ + 1) 
for all t < t2 — (J , a > 0. 
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Proof. Note that the function C^(^) in (14. 9p is C°° in (—00, t2~cr), and that each 
Uk, k = 0,1, . . . , given in fl3.27p . is continuous at t = 0. Hence the expansions 
in (13. 4 p and (13. 5 P hold uniformly with respect to t in any compact subinterval of 
(—00, ^2); and the coefficients Gs{C) and Hs{() in these expansions are bounded 
for all finite t away from ^2- Thus, to prove the two estimates in (15. ip . it suffices 
to show that the following two limits exist: 

(5.3) lim li,„ m 

t— > — 00 t t— — 00 t 

Replacing Zj, by and and choosing positive sign in (13. 2p . we obtain 

(5.4) J, {{N + l)C(t+)} = MNOg{c, ^) + J.+i{NOh(c, I 
and 

(5.5) n {{N + l)aU)} = Y,{N()g(c, 1) + Y^^,{N()h((, i) . 

Upon solving the last two equations for G{(, -^) and H{(, j^), we have 
(5.6) 

g{c, ^) =^{Ju^iiNOY. [{N + l)C(t+)] - n+i(iVC) J. [{N + l)C(t+)] } 

and 
(5.7) 



h(c,- 



{Y^iNOJ. [{N + l)C(t+)] - MNQY, [{N + l)C(t+)] }, 



2 

where we have made use of the identity 

2 

J^+i{x)Y^{x) - J^{x)Y^+i{x) = — ; 

TCX 

see [15], p. 244] and [16, eq. (10.5.2)]. Recall the relations between the Bessel 
functions and modified Bessel functions fibl pp. 60 & 251] 

(5.8) Mix) = e'^'I^Ux), Y,{ix) = ^x) - -e-'^^^K^ix), 

IT 

and let 

(5.9) ^ = -iNCit), e' = -i(iV+l)C(t+). 
Then, (15. 6p can be written as 

(5.10) G^cj^) = ^[hiOK.MO + KiOlu+iiO]- 

Recall the asymptotic expansions of the modified Bessel functions for large 
argument 

(5.11) ^"W-p^E(-i)f 
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and 



(5.12) 



KJx) 



/ TT \ 1/2 



s=0 



where /xq = 1; see ^ pp. 250-251, 238]. A combination of (I5:T0|) - (I5:T2| ) gives 



(5.13) 



2 If 



s=0 



1=0 



~ oo / ^ \ I 



+e«'EM)f Ef I 

s=0 ^ 1=0 ^ 



where JIq = 1 and 

Jls{iy) = fisii" + 1), s > 1. 
It follows from IK9\} and ( liTTOl) that as t ^ -oo 



C^\t) = O {\t\^^) , e>2, 



(5.14) 
and 

(5.15) c^^(t) = c (log^^ , o<e<2oTe<o. 

Using (14.91) . f l4.10p and fl5.14p . a straightforward calculation gives 



<t + + v/(«ot + Pof 



(5.16) 



+ 



2 

[s-l)(l 



iV Jo v/(«^t(l + ^)-^ + /3^)2-4 



rfs := ^ 



where 6o{t) = log[(a[)t + + yKtT^^F^)/2] and 



(5.17) 



6,{t) 



Oj {\t\'^) as t -> -oo. 



J(J + 1) 

Note: in (15.171) we have used a subscript "j" to indicate that each C-term 
depends on j, j = 1,2, . . . . Hence, 



(5.1^ 



For each j > 1, the coefficient of is (log ^\t\) for large negative t. We 
define 



exp 



exp 



oo 



j(j + l) 



j=0 



aj(t) 



-E 
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and 



for each 



It is readily verified that ao(^) = ^o(^) = Q,o(^) = 1, cijit) = ^jilA 
j > I and each / > 0, and other bj{t), cij(t) are also bounded for t < t2 — c. 
With these notations we have from f l5.13p that 



2G(C,^) ~exp 



aj it) 



X 



s=0 



+ exp 



«=0 



AT' 



m=0 



{-lye 1 



j=0 



C/+im(^) 



X 



E(-i)':^('r'w)'-E^(T'w)'E m- 



s=0 



m=0 



which in turn gives 



-(N+i)e 



p=0 



NP 



-(N+l)9 



where Wifi(t) = W2fi(t) = 1 and for p > 1 

j+s+l+rn=p 
j+s+l+m=p 

Note that for each p > 1, we have wi^p = O (1) and W2,p = O (1). If we set 



then we further obtain 



(5.19) 



-(Af+i)e oo 



E 

p=0 



ATP ' 
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where ro(t) = 1 and Tp(t) = O (1) for each p > 1. Let 

-(7V+1)6» 



(5.20) 



exp 



^-0(iV + l)log( 1 + - 



A combination of f l5.19p and fl5.20p yields 
(5.21) G{(,^)~'^±l-±...Mt), 



s=0 



j=0 



and a comparison of (13. 4p and (I5.2ip gives 
(5.22) 



G.(C) = ^^^[^^ + 0(l) 



The existence of the first hmit in (15. 3p is thus proved. The existence of the 
second hmit in (15. Sp can be shown in a similar manner, and we have thus 
established the two estimates in (15. ip . The proofs of the two inequalities in 
(15.20 are very similar to the ones given above, and hence will not be repeated 
here. This completes the proof of Lemma [5. II □ 



Using the estimates proved in Lemma 15.11 we can provide some bounds 
for the coefficient functions As{() and -Bs(C) in (14. ip . We state the result as 
follows. 



Lemma 5.2. For s, j 

such that 



0,1,2,... and t < 0, there exists a constant Msj 
t^'-^i-Ho/Ct-D^MC)] < Msjil + |trf ) 



(5.23) 
and 

(5.24) \t^+^{-Ho/Ct'D^B,{0\ < M,,,(l + \t\~^), 
where denotes the j-th derivative with respect to t. 

Proof. Our approach is based on mathematical induction. Since Gs{C)y 
Hs{C), Ls{C) and Ks{C) are all analytic functions in the left half plane Ret < 0, 
and since the estimates in Lemma ISTT] also hold for complex t in this half plane, 
an application of Cauchy's integral formula shows that for s, j = 0, 1, 2, . . . and 
t < 0, there exist constants Nj and Nsj such that 

(5.25) \tW^(:\ < iV,(l + Id), 

(5.26) \tW^Gs{0\ <Ns,ja + \t\), 

(5.27) \fD^H,{C)\ <Nsjil + \t\), 

(5.28) \tW^Ks{0\ <iV,,,(l + |t|), 

(5.29) \t'D^L,{0\ <iV,,,(l + |t|). 
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When s = 0, the estimates in f l5.23p and (15 .24^ follow from f l4.25p for j = 0. 
This, together with fHTTD and dHHD, implies that fl^:^ and hold for 

j = 1. Furthermore, by differentiating f l4.17p and (14.181) j times with respect 
to t, one can show that (I5.23P and (I5.24p also hold for j > 2. Now assume 
these estimates hold for s = 0,1, ... ,p — 2, p > 2. We shall show that they 
are valid for s = p — 1. Using (I5.25p - (l5.29p . it can be shown from (I4.15P and 
(|4J6|) that 



\t^{-Ho/Cr'fp-iit)\ < Cp,o(l + |t|)(l + \t\ 



p-2 . 





p-2 . 

e 



(5.30) 
and 

(5.31) |t^(-iJo/C)^^?p-iWI < Cp,o(l + |t|)(l + |t| 

where Cp^ is a positive constant. Combining the last two inequalities with 

KTHf leads to 

(5.32) 



p-i 



C 



p-i 

S e 



1 OHoiCis)) 



A{s)fp^i{s)ds 



<\t\ 



C + CpX) 



, p-1- 



p-1 _ 



.p-2 



)ds 



[T 



+ re 



<Mp_i,o(l 



.P-i , 



Here we have made use of (14.191) . In a similar manner, we have 



(5.33) 



\t^{-Ho/Cr^Bp.,iC)\ < Mp„i,o(l + |t| 



_p-i , 



Hence, and hold for s = p - 1 and j = 0. Using (HTTD and 

(I4.18p . one can also show that (I5.23P and (15.241) are valid for s = p — 1 and 
j = 1, 2, . . . . This completes the proof of Lemma [5.21 by induction. □ 

6. Asymptotic nature of the expansion 

Since x is fixed in the recurrence relation (11.71) . we let Wu{x) := Y^{x)—iJ,^{x) 
and choose 

Z^{NC) = x^J^{NC) and Z^{NC) = x^W^{NC) 
in (14. ip to yield two hnearly independent solutions 



(6.1) 



AC 



4 - {a'ot + P'^ 



'\2 



N2 



s=0 



MO 



s=0 
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and 



(6.2) 



s=0 



s=0 



where ^(C) = m [- HoiQ / C]'^ MQ , 5.(0 = [-//o(C)/C]^i?.(C), and 

(6.3) ^o(C) = l, 5o(C)=0. 

Here we have made use of ( KWf . (14241) and ( 105ll . From K9\\ and ( HTTOj) . one 
can get the behavior of C(t) as t — > -co; cf. f l5.14p and f IS.lSp . With this, it 
follows from Cauchy's integral formula that there is a constant Ck such that 

\t''D''A{t)\<Ck\A{t)\, t<0. 

Applying Leibniz's rule to the products As(C) = ^(^)^s(C) aiid -Bs(C) = 
A{t)Bs{C), the estimates in (lOSj) and (ISlMj) give 

(6.4) |t^Da,(C)| <M.,,(l + |tr^) 
and 

(6.5) \tW^B,iC)\<Msjil + \t\-i), 

where Ms,j is some other constant. It readily follows from fl6.3p that Aq{() 

and -Bo(C) are bounded for t < t2 — a. Assume that As{() and -Bs(C) are 
bounded for t < t2 — c and s < p. Using an argument similar to that given 
for fICTD - flOSD . we have from ( KT5\f and ( KWf that 

|t^(-iJo/C)^/p+i(t)l <Cp+i,o(l + |t|), 

\t^{-Ho/Ohp+i{t)\ < Cp+i,o(l+ |t|). 
As a consequence, it can be shown that fl5.23p and f l5.24p can be improved to 
read 

and 

for t < and s,j > 0. In particular, these two estimates hold for j = 
and s = p + 1] that is, Ap+i{Q and Bp^i{() are bounded for t < 0. By the 
continuity of C(^), Gsi() and Hs{C) we conclude that y4p+i(0 and Bp^i{() are 
also bounded for t <t2 — a, thus completing the induction argument. 

Now we shall show that (11. 7p has two linearly independent solutions Pn{x) 
and Qnix) with (16. ip and (16.20 being their uniform asymptotic expansions for 
t<t2-cT. We need to show that fl^:^ and fl^:^ hold for t <t2-(T. 
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For convenience, we introduce the notations 

4C' \K,i^MO 



s=0 



(6.8) := MNC)A,[C ^) + J.+i(iVC)5,(c, ^ 
and 

(6.9) <(x) := W^,(iVC)A,(c,^) +W^.+i(iVC)i?,(c,^). 
By Lemma [3.11 we have 

(6.10) Ci(x) - {A^x + B^Kix) + r^,(x) - ^'^""^ 



NP+I ' 



where x = N^t and the inhomogeneous term is given by 
(6.11) ^ = J.(iVC)Fi,„(x) + J.+i(iVC)F2,„(a;) 

with 
(6.12) 

F.A.) =g(c i)^,(c(tj.^) +G(c.-^)^,(C(i-).^ 

+ <C. ^)b,(c(*.). ^) - i(c, 4)i',(c(«-). ^ 

-*(*-^)-4'(<-i) 

and 
(6,13) 

F.,.(.) ='^-(c4)«''(';(*+)-FTl) +'^'('^'4)^''(«'-''7v^) 

-*('-^)^^(<-^)^ 

cf. (gSD and fOj) . Note that the series 



s=0 s=0 
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in f l4.2p are formal solutions of ( I4.5p and f l4.6p . Since Ap (^(, j^) and Bp j^) 
can be rewritten as 



(6-14) ^p(C,^):=x^$:^ and Bp{c, j^) := ^1 

s=0 s=0 

With A:(C) = MO, B:{0 = 5.(0, for ^ < p and A^Q = 0, i?:(C) = for s > 
p, terms with powers of less than or equal to p+1 in the expansions 
and F2^n{x) all vanish, i.e., if we write {N^t)~^ Fi^n{x) ~ Xl^o ^'^'^ 
(iV^t)-WF2,„(x) ~ YZo f2,s{t)N-\ then /i,,(t)=/2,.(t) = for < s < p + 1. 
(For s = 0, 1, . . . ,p, all we need are recurrence relations obtained from fl4.5p 
and (14.61) : for s = p + 1, there are extra terms and we have to show that these 
terms indeed all cancel out.) Using (I5.6l) - (l5.12p . it can be shown that there is 
a constant C such that 



g(c, ^) I < c{i + \t\), \h{c, ^) I < C{1 + \t\) 



and 



(c, ^) I < c(i + \t\), |l(c, ^) I < c(i + |t|) 

for all t < t2 — (J. Recall that As{() and -Bs(C) are bounded in t < t2 — cr for 
s < p. Hence by Lemma 15. it can be proved that there exists a constant Cp 
such that 



(6.15) 
and 
(6.16) 



\t^i-H,/C)HN't)-^F,^^iNh)\ < Cpil + |t|)/iVP+2 



for all t < t2 — cr. Then, the last two inequalities, together with (16. lip , lead to 

AC 



(6.17) mNH)\ < Cpil + \t\) 



[\MNC)\ + \J.MNC)\] 



4 - ia'.t + 

for t <t2 — CF and for some positive constant Cp. Similarly, we have 



(6.18) 
with 



Sn+li.x) - iAnX + Bn)slix) + S 



n—l\ 



(6.19) \Sl{Nh)\<Cpil + \t\) 



AC 



4 



[|i^.(ivc)| + |w^.+i(ivc)|]. 



We now establish the existence of two solutions Pnix) and Qnix) of ( II. 7p 
satisfying 

(6.20) P„(x) ~ r°(x) and Q„(x) ~ s°(x) 
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as n — )■ oo for x = N^t with any fixed t < t2 — u. In view of (15. Sp . by using 
flS.lip and (15.121) it is easily verified that 

and 



as n ^ oo for any fixed t < 0. From ([321, (I32ZD, (E^H]) and (gS]), it follows 
that 



(iV + l)C(t+)-iVC(t)~-log 

1 

as n — oo. If Pri(x) and (5„(x) are two solutions of (I1.7P satisfying (I6.20p . 
then we have from (lOOD - flO^ 

(6.23) Pn+l{x)Qn{x) - Pn{x)Qn+l{x) = Pn+2{x)Qn+l{x) - Pn+l{x)Qn+2{x) 

and 

Pn+lix)Qnix) - Pn{x)Qn+lix) 

(6.24) 2 

= lim [^!^i+i(a;)s[|,(x) - r[^(a;)s[^+i(x)] = - 

m— ^oo 71" 

for fixed t < 0. Examining the behavior of the Bessel functions, one can show 
that (I6.24P also holds for fixed t > 0. This, in particular, demonstrates that 
Pn{x) and Qn{x) are two linearly independent solutions. 
Now define 

(6.25) el{x) := P„(x) - <(x) and 5P{x) := Q„(x) - s^(x). 

We first show that the existence of Qn{x) to (II. 7p satisfying (I6.20p is equivalent 
to the existence of S^{x) to the summation formula 

APr ^ _ K(a;)'^i(a;) -rP(x)s^(x)]5^(a;) 



i=n+l [Sn+l{X)rn[X) - Sn[x)r^^^{x)]{j + TqY 

6.26) 

^ ' K{x)R^^{x)-r^^{x)S^{xW^{x) 



3 



j=n+l [■'n+l(a;)'"n(x) - (x)r^+i (x)] (j + To)P+2 



Prom (ll.7p and (16.180 . we obtain 

(6.27) 5^^,(x) - {A^x + B^)6l{x) + 6l_,{x) = 
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Coupling f l6.10p and f l6.27p gives 

(6-28) g„^i(x)r^+i(x) + /?::^,(3:)Ci(^) 

(A^ + l)P+3/2 

which leads to 

(6.29) S^{x)r^{x) + 

If (5„(x) satisfies fl6.20p . then S^{x) = o(s°(x)) as n — )■ oo. Since = 
(5^(a;) + — s^(a;) by fl6.25p . as well as As{() and -Bs(C) are bounded for 

s > 0, on account of fl6.8p we also have S^{x) = o(s°(x)) as n — )• oo. In view 
of (ICT]) and we have 

(6.30) C+i(a;)5^+2 W - '^m+2(a;)C+i(a;) ^ 
as m — 7- oo. Hence, by letting m — ?► oo in (16. 29 p . we obtain 

(6.31) r^(x)5^+i(x) -<+i(a;)5^(x) = 2^ ^ ^ • 

In a similar manner, we also have 

(6.32) <(x)Ci(-) - sl,,{x)5l{x) = ± ^^^"^f f ^^.t/f ^'^'^ 

j=n+l 

and 

(6.33) r:,,{x)s:{x)-sl,,{xK{x) = l+ ^^^"^^.^"^^^^^^^'^^"^ 

Then (I6.26P follows from (I6.3ip and fl6.32p . The existence of a solution {S^{x)}'^^^ 
to fl6.26p is proved by using the method of successive approximation. Starting 
with S^q{x) = 0, we define S^i^{x) by 



j=n+i K+i{x)rZ{x) - s^(x)r^+i(x)](j + ro)P+2 



:^)d(x)](j + ro)^+^ 



for k > 1 recursively. We shall show that for fixed p and sufficiently large but 
also fixed n, the sequence {S^ k{x)}k>o is convergent as A; —t- 00. Since As{() 
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and -Bs(C) are bounded for t < t2 — a, it follows from f l6.6p -( l6l9l) that 

4^2 



(6.35) 
and 
(6.36) 



\sUx)\ < cm 



4 - {a',t + /3^)2 
AC 



[\UNQ\ + \J,^,{NQ\] 



[|H^.(iVC)| + |l^,+i(iVC)|] 



for some positive constant C. Furthermore, by virtue of the behaviors of Jy 
and Wy, we have from (ICTTl) . ( I619|) . (16351) and (lOeD that 



(6.37) 



|i?^(x)sP(x)| < M'ATs and |^p(xX(x)| < M' n'^ 



for some positive constant M' . It then follows from f l6.33p that 



sl+i{x)rl{x) - r^+i(xX(x) + 



TT 



< 



2M' 1 



p N'P 



Since the right-hand side tends to zero, this estimate gives 

(6.38) \sl^,{x)rl{x)-rl_,,{x)sl{x)\>- 

7r 

for large n. A combination of f l6.34p and fl6.35p -f lB.38p yields 

[|.^„(x)<(a;)| + |r^(x)4(x)|]|5^(x)| 



j=n+l 



<M'' 



4 - {a',t + %y 



(j + ro)^+i 



j=n+l 



A^2 



(j + ro)P+ 



— [|w-.(ivc)| + |iy.+i(ivc)|]: 



where we have also made use of the monotonicity properties and the asymptotic 
behaviors of the modified Bessel functions and Ky] see fl5.8p . f l5.1ip and 
flHTT^ . Hence, 



M" 1 



4C^ 



(6.39) |<i(x)|< 

Similarly, we can prove by induction that 



[\Wy{NC)\ + \Wy+,{NC)\]. 



(6.40) 



< 



M" 



from which it follows that 



4C^ 



(6.41) 



k 



[|w^.(ivc)| + |i^.+i(ivc)|]. 



m=l 
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converges as k ^ oo, for all n > 2M" — tq. Clearly, the limiting function S^{x) 
satisfies (16.26 p . Thus Qn{x) = s^(x) + 5^ (a;) is a solution of (II. 7p satisfying 
(I6.20p . Furthermore, we have from (I6.39p - (l6.4ip that 

Recall that As{() and -Bs(C) are bounded for t < t2 — u- By taking an extra 
term in the expansion in (16. 2p . we have 



(6.42) X [\W,{NO\\A,+,{0\ + |W^,+i(iVC)||5,+i(C)|] 



< 



[|w^.(ivc)| + |w^,+i(ivc)|] 



for some positive constants and Mp, and (I2.2ip follows. (Note that there is 

an extra factor of A^^ outside the square brackets in (I2.19p . so the error terms 
in (I2.2ip and (I6.42P have a slightly different meaning.) 
In a similar manner as we have done for S^{x), it can be proved that 

(6.43) < ^ [t^Stw) ' ^ 

and hence (I2.20p follows. This completes the proof of Theorem [TJ 



7. An EXAMPLE 

As an illustration, we consider the orthogonal polynomials Pn{x) associated 
with the Laguerre-type weights w{x) = exp{—qjnX"^), x > 0, a > — 1, 
Qm > 0. This weight is a Freud weight on the half line; the asymptotics of 
Pn{x) have been investigated by Vanlessen [19] via a Riemann-Hilbert method. 
For the Freud weight on the whole real line, the asymptotics of the associated 
polynomials can be found in [9l [HI |26]. The notation Pn{x) = jnx'^ + • ■ ■ with 
leading coefficient 7„ > denotes the n-th degree orthonormal polynomial 
with respect to w{x). It is known that the sequence of Pn{x) satisfies the 
three-term recurrence relation 

(7.1) xpn{x) = bnPn+i{x) + a„p„(a;) + 6„-iP„_i(a;), n = 1, 2, . . . , 

with p^i{x) = and ^0(2^) = 7o > 0. The recurrence coefficients a„ and bn-i 
have the asymptotic expansions 

(T.2) 6„_.~n^rJUfi + f| + 0(i)| 
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and 

, 1 f 1 a + 1 1 (ia / 1 

where and d2 are some constants and 

(7.4) Tm = -K^qm JJ 



see [19]. We define a sequence {-ft'n} by 

n + 21 \ 2k 6„+2/+i 



llLVn + 2/ + 2/ 



n + 2/ + 2y h^+2i J 
Note that -ft'„+i/-ft'„„i = hn-i/hn and ~ n~2^. We then put 

(7.5) = = P„(x) ^ (-l)"Kx)]V„(x)^ 
The three-term recurrence relation in fl7.ll) now becomes 

(7.6) Vn+l{x) + Vn-l{x) = {A^X + Bn)Vn{x). 

A slight computation gives the following asymptotic expansions 
Kn ,,11 4m + 2am - 1 1 , f I 



(7-7) -W^-^+o ^2 2 + ^ - 

2m n 8m^ n"^ \ n'^ 



(7.8) ~ n-^r-^ {-4 + ^±^1 + o(\ 

and 

[Aw? — Am + 1)0^ — rn? 1 
7.9 i?n~2 + ^ -^^ ^ + 3 

In terms of our notations, = ^ and tq = ^j^. Set = n + tq. We then have 



00 



s=0 s=0 

where 

1 / ^/ ^/ ^/ (2m — l)^a^ — m^ 

«[, = -4r-\ a; = 0, /3^ = 2, /^J = 0, ^ ^^^^^ . 

Let X = N^t. To simplify matters, we further set t = and u;„ = N^rm, so 
that X = cOnZ. Since ^ = 7^ 0, 2 and /Sg = 2, the two transition points are 
ti = and ^2 = ■'"m! cf. fll.llj) . We can now apply Theorem [H In terms of the 
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new scaled variable 2;, the two transition points are = and Z2 = 1- From 

f l2.15p . it can be shown that u = a, and by (14. 9 p 

(7.11) 

C(^) = arccos(l - 2z) + — ^2Fx ( 1, 1 



2m 



z < 1 



a. 



Remark 2. For any given integer m > 0, the hypergeometric function in 

(7.11) is a polynomial of degree (m — 1). For example, when m = 1, it is 1; 
when m = 2, it is 2z + 1; and when m = 3, it is + |z + 1. 

From Theorem [H it follows that we have two linearly independent solutions 
Pn{x) and Qn{x) given in fl2A8|l and fl219D . Hence, 

(7.12) VnicOnZ) = Ci{x)Pn{x) + C2{x)Qn{x), 

where Ci{x) and C2{x) are functions depending only on x. To determine Ci{x) 
and C2{x), we note that in terms of our notations, equation (2.12) in ^19j reads 

'2r[{2z-l)+2^z{z-l)X 



(7.13) [w{(5nZ)]-^Pn{(3nZ) 



2 



2z^{z 



1)^ 



exp{-inC(2;)}, 



where /5„ = n^rm- This asymptotic formula holds uniformly for z in any 
complex subset K C C\ [0, 1]. Note that our uniform asymptotic expansion of 
Vn{ujnz) f l7.12p is also valid in a neighbourhood of any subinterval of (—00, 1 — 
a] in the complex 2;-plane. Since = (1 + ro/n)m/3„, matching these two 
asymptotic formulas in an overlapping region, and using the behaviors of J,y{x) 
and Yi,{x) for large x, we obtain Ci{x) = 1 and C2{x) = 0. That is 



(7.14) 



^(1-^1 



s=0 



MO 



s=0 



where Ao{C) = 1, 5o(C) = 0, and 

M 



(7.15) 



NP+ 



'^[\UNO\ + \Ja-,l{NO\] 



for all z < 1 — a. We can also use the result in [21] to derive a uniform 
asymptotic expansion for Vn{x) near the turning point Z2 = 1, and we have 



i-irVnioOnZ) = iV6 



(7.16) 



z{z - 1) 
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for z > a > 0, where ^0(77) = 1, Bq{vi) = 0, 



(7.17) ^-r]{z)Y = arccos(2^ - 1) - \^ _ ^ ■ 2^1 f 1, 1 - m; - - m; z 

for 0" < 2; < 1, and ri{z) analytically continued to C \ {(— oo,0] U [1, — ioo)} 
such that ri{z) > foT z > 1. The error estimation is given by 

(7-18) 1^1 < ]^Ai(ivl.^) 

with the modulus function Ai(A^i?7) as defined in [21[ eq. (7.10)]. Note that 
the two results in f l7.14p and fl7.16p together cover the whole real line. 

Remark 3. In the special case of m = 1 and qm = 1, we get two uniform 
asymptotic expansions for Laguerre polynomials, which agree with the results 
obtained by steepest decent method for integrals |TT] or WKB approximations 
for differential equations [151 Chaps. 11 and 12]. 

For orthogonal polynomials associated with the weight exp(— (5(x)), x > 
0, a > — 1 and Q{x) is a polynomial of m-th degree with positive leading 
coefficient, the asymptotic expansions of the recurrence coefficients will be in 
powers of l/n™, instead of in powers series of 1/n; cf. M.^) . In such a 
case, one can modify the method provided in this paper to get a pair of linearly 
independent solutions to the three-term recurrence relation given in TheoremUl 
which are also in terms of Bessel functions near the transition point ti = 0. 
However, the asymptotic expansions for these solutions are on longer in powers 
of 1/N, but in powers ofl/N^, This result will lead to Bessel-type asymptotic 
expansions given by Vanlessen in 



As another example, confluent hypergeometric functions M(a, b; x) and r(a— 
b + 1) ■ U{a,b; x) satisfy the three-term recurrence relation in a: 

(7.19) ay{a + l,b;x) + {a — b)y{a — 1, b; x) = (2a — 6 + x)y{a, b; x). 

A straightforward calculation shows that equation fl7.19p corresponds to our 
case for ^ = 1, ao = 1 and Pq = 2. A direct application of Theorem [1] yields 
Bessel-type expansions as a — )■ +00. Replacing a by —a in fl7.19p . one can 
obtain again, by the main theorem, Bessel-type expansions for M(a, b; x) or 
U (a, b; x) as a — )■ —00; see [121 §13.8(iii)]. 
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